Extraction Of Ag/g From Hermes Data On 
Inclusive Charged Hadrons 
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Abstract. Hermes has measured longitudinal double spin asymmetries as a function of transverse 
momentum pj using charged inclusive hadrons from electroproduction off a deuterium target. At 
Pt > 1 Ge V, the asymmetries are sensitive to the spin dependent gluon distribution Ag. To extract the 
gluon polarization Ag/g, information on the background asymmetry and the subprocess kinematics 
has been obtained from a Leading Order Monte Carlo model. Values for Ag/g have been calculated 
both as a function of the measured pj and x, using two different methods, in the region pj > 1.05 
GeV. 
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INTRODUCTION 



A direct, model dependent extraction of Ag/g has been performed by Hermes (H], SMC 
[0] and Compass [0, 13] for different channels and data sets. This report presents a 
refined extraction method, using the high statistics data sample of antitagged (vetoed 
by electrons or positrons), inclusive charged hadrons. A detailed study was performed 
to estimate the systematic error arising within the model which uses Pythia 6.2 [5], and 
parametrizations of spin dependent parton distributions of the nucleon and of the photon. 



EXPERIMENTAL DATA 



^ ■ The data sample used for this analysis was collected using the Hermes spectrometer. 

Charged inclusive hadrons were selected from events where neither a positron nor an 
electron were detected. The asymmetries were calculated as 

where are the number of hadrons detected with beam and target spins parallel (an- 

tiparallel), are the corresponding integrated luminosities and P~K~~) the integrated 

luminosities weighted with the product of beam and target polarizations. The transverse 
momentum was calculated with respect to the beam axis. The results for positive and 
negative hadrons from proton and deuteron targets are shown in figured] The asymme- 
tries have not been corrected for acceptance and trigger efficiency. It was confirmed that 
the trigger efficiency does neither introduce a significant bias to the asymmetries nor to 
the final results. Also shown in the figure are the asymmetries expected from the model, 
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FIGURE 1. Measured asymmetries for antitagged inclusive charged hadrons. The top row shows the 
asymmetries for proton- , the bottom row for deuteron target. The curves show the asymmetries calculated 
using Monte Carlo and spin dependent quark distributions and the assumptions Ag/g(x) — —1,0 and +1 
(from top to bottom). 



using the assumptions Ag/g(x) = —1,0 and +1 (lines from top to bottom) over the full 
covered x-range. More detailed information on how the model asymmetries were calcu- 
lated will be given below. The differences between the measured asymmetries are due to 
the quarks, and are well described by the curves at low pr, where the contribution from 
gluons is negligible. 



INTERPRETATION IN TERMS OF THE GLUON POLARIZATION 

The measured asymmetries arise from a superposition of different subprocesses con- 
tributing to the production of hadrons at a given measured pj. In order to decompose 
the asymmetries and extract the signal from processes initiated by a hard gluon, the 
asymmetries and relative contributions of the background processes have to be known 
as well as the hard subprocess kinematics of the signal processes. This information was 
obtained from a simulation of the data using the Pythia 6.2 Monte Carlo program and a 
model of the Hermes detector. 



Data - Monte Carlo Comparison 

The Vector Meson Dominance (VMD) Model in Pythia was adapted to reproduce the 
observed exclusive p° cross section [|6j]. The fragmentation process simulated in Jetset 
was tuned to multiplicities of identified hadrons at Q 2 > 1 GeV 2 measured at Hermes 
10]. For Q 2 > 0. 1 GeV 2 , the observed semiinclusive cross sections agree typically within 
15% in variables integrated over pj. In contrast to this generally good agreement, 
the observed cross sections do not agree vs. pj. For pj > 0.7 GeV, the Monte Carlo 



underestimates the data by factors 2 to 4. The disagreement is most likely due to large 
NLO corrections as calculated in Ref. [8]. The LO result for Ag/g might therefore also 
be subject to large NLO corrections. 



Subprocess Fractions 

The soft background processes from exclusive and diffractive VMD as well as non- 
diffractive VMD ("low-py") are suppressed at high pj, although the "low-p^" process 
still contributes significantly for pj < 1-5 GeV. The quark initiated hard QCD processes 
together contribute less than 20% at pj > 1 GeV. At pj > 1.5 GeV, the LO DIS process 
dominates the cross section. Hadrons from this process originate predominantly from 
events with a large lepton scattering angle, where the transverse momentum calculated 
with respect to the beam axis overestimates the transverse momentum in the center of 
mass frame. The signal processes are Photon-Gluon-Fusion (PGF) and the gluon initi- 
ated 2 — > 2 (resolved photon) processes, each contributing 10-20% in the relevant p T 
range. 



EXTRACTION METHODS AND RESULTS 

The background asymmetries for hard subprocesses were estimated using the MC infor- 
mation on the particle types and subprocess kinematics with the nucleon PDFs from [@] 



and the photon PDFs from HOD, where the average of the maximal and minimal scenarios 
was taken. The asymmetry for exclusive VMD was set to 0, that of the "low-pr" pro- 
cess was set to A low ~ PT = g\/F\, using an extrapolation of a fit to world data to lower 
x ~ 10~ 4 . By subtracting the background asymmetry weighted with the background 
fraction from the measured asymmetry, A™ eas — RbgA bg = R s i g A slg the signal asymme- 
try can be obtained. The signal asymmetry contains a convolution of Ag(x)/g(x) with 
the (polarized) hard subprocess cross section over the jc-range covered by the data. Two 
methods have been applied to extract the average (Ag / g) (pj) from this asymmetry using 
different assumptions on the shape of Ag(x)/g(x). Method I assumes that Ag(x)/g(x) is 
essentially constant in the relevant jc-range. Then, (Ag/g)(/?r) can be found by solving 
the equation 

Aj""-^^^)^) (2) 

for each bin in pj. Here, a is the hard subprocess asymmetry, and Ap / p is the 
polarization of partons in the resolved photon. 

In Method II, a functional form is assumed for Ag(x)/g(x) which is used to calcu- 
late the integral A slg for each pj bin. The functional parameter(s) are determined by 
minimizing % 2 for the difference A™ eas — RbgA bg — R S i g A s ' 8 using all bins in pj. 

Values for (Ag/g)(pr) have been obtained from the deuterium data on charge com- 
bined hadrons, for 4 bins in pj between 1.05 and 2.5 GeV. The results are shown in 
figure [2] for both Methods. Note that, in Method II, the results and errors are corre- 
lated through the fit function. Method I was used to confirm the overall consistency 
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FIGURE 2. Results for (Ag/g)(p T ) from Method I (left) and Method II (right) for the anti tagged 
deuterium data. The inner error bars (band) correspond to the statistical errors, the outer error bars (band) 
to the total errors. 



between different independent data sets from proton and deuteron targets and positive 
and negative hadrons. The experimental systematic error is approximately 14% and 
arises from the uncertainties in the beam and target polarization measurements. It is 
small compared to the model uncertainty which was estimated by varying the Pythia 
model parameters, the unpolarized PDFs in the MC generation, the polarized PDFs in 
the asymmetry calculation and the assumption used for the asymmetry of the "low-pr" 
process. For Method II, the functional form was also varied. No error was assigned 
on the Pythia model itself, and, because this is a leading order approach, also no error 
was assigned to account for NLO corrections. Method II allows to determine the av- 
erage x of the measurement, and by integrating over 1.05 < pr < 2.5 GeV a value of 
Ag/g = 0^1 ±0m4{stat) ±0M0(sys - exp)t° \l 7 5 {sys - Models) has been obtained 
at (x) = 0.22 and (ji 2 ) = 1.35 GeV 2 . 
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